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Inelastic neutron scattering is used to study the finite-temperature scaling behavior of the local 
dynamic structure factor in the quasi-one-dimensional quantum antiferromagnet NTENP (Ni(A', N'- 
bis(3-aminopropyl)propane-l,3-diamine)(/r-N02)C104), at its field-induced Ising quantum critical 
point. The validity and the limitations of the theoretically predicted scaling relations are tested. 


I. INTRODUCTION 

The Ising chain in a transverse field is one of the 
most-studied model systems for quantum phase transi¬ 
tions thanks to its simplicity and the possibility to ob¬ 
tain exact results.lill^It features a quantum critical point 
(QCP) at a certain value of magnetic field He where the 
ground state changes from a quasi-classical ordered to 
a disordered quantum-paramagnetic phase.^^ Early the¬ 
oretical studies examined the ground-state energy and 
spin correlations at zero temperature, T = O.^i^^More re¬ 
cent work was able to derive explicit expressions for the 
finite-tempe ratu re correlation functions precisely at the 
critical field,!^ and even in the qua ntum -paramagnetic 
and the magnetically ordered phases.^^ Independent of 
the details of the actual Hamiltonian, in the quantum 
critical regime, the spin correlations follow a universal 
scaling law. Such universality is particularly interesting 
for experimentalists for two reasons. On one hand, all 
predictions are expected to hold for any material featur¬ 
ing a QCP in the universality class of Ising chains in a 
transverse field, irrespective of any material-dependent 
quantities. On the other hand, spin correlations can 
be directly and accurately probed by neutron scattering 
experiments enabling an immediate comparison between 
experiment and theory. 

Unfortunately, only few experimental realizations of 
spin-chain materials with strong uniaxial anisotropy and 
accessible field-induced QCPs are known. Furthermore, 
most exhibit strong interchain coupling, as in the case of 
CoNbaOfplor CaaCoaOfpl. Since strong interchain cou¬ 
pling leads to long-range order at finite temperatures al¬ 
ready below He, these materials are only of limited use 
for studying the QCP. Fortunately, the one-dimensional 
(ID) quantum Ising universality class - which is equiv¬ 
alent to the two-dimensional (2D) thermodynamic Ising 
universality class - is also realized in the field-induced 
transition of Heisenberg S' = 1 chains with easy-plane 


anisotropy (Haniso = D > 0) and with 

a m agn etic field applied in an axial-symmetry-breaking 
integer spin is necessary in combination 
with an easy-plane anisotropy in order to map the system 
to the one-dimensional quantum Ising model. For this 
class of Hamiltonians, excellent experimental realizations 
are known, and have been studied in great detail. The 
maybe best-known examples are Ni-based linear- chain 
azide complexes, e.g. NENPG^, NMOAfC^l^ NDMAPfSES 
or NDMA^. 

The objective of the present work is to exper¬ 
imentally test the finite-temperature scaling of the 
spin correlation function at the ID Ising QCP 
with neutron scattering. The target compound of 
this study is the antiferromagnetic (AF) S' = I 
chain material Ni(fV, fV'-bis(3-aminopropyl)propane-l,3- 
diamine)(/r-N 02 )C 104 (NTENP). The compound is well- 
characterized. At room temperature, NTENP is a tri¬ 
clinic system with space group PI, lattice constants 
[a,b,c] = [10.747(3), 9.413(2), 8.789(2)] A and angles 
[a,/3,7] = [95.52(2)°, 108.98(3)°, 106.83(3)°].E3 It fea¬ 
tures bond-alternating AF exchange with exchange con¬ 
stants Ji = 2.1 meV and J 2 = 4.7 meV as well as single¬ 
ion anisotropy, P = 1.2 meV.^^SHU] Interchain interac¬ 
tion are only weak, being of the order of 0.005Ji.l^ The 
magnetic properties of NTENP can be described by the 
following leading terms in the Hamiltonian 

N/2 N 

n = Y, {JiS 2 i-iS 2 i + J2S2.S2.+1) + Y,D (sf)". ( 1 ) 

i—1 i—1 

The chains along the crystallographic a axis are formed 
by Ni^+ ions which interact via NO^ groups.^ The crys¬ 
tal structure is shown in Fig.[2 The material has a spin- 
dimerized non-magnetic ground state and the excitation 
spectrum is, accordingly, gapped. No magnetic order can 
be observed down to 50 mK in the absence of a mag¬ 
netic field. The first excited states are an 5 = 1 triplet, 
split by single-ion anisotropy into a lower-lying doublet 
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with a ga p of 1.06 meV and a higher-energy singlet at 

1.96 meVApplying an external magnetic field de¬ 
creases the energy gap in one member of the doublet due 
to the Zeeman effect. For an axial-symmetry-preserving 
magnetic field applied along the hard axis, i.e. perpen¬ 
dicular to the easy plane, a transition to the gapless 
Tomonaga-Luttinger spin liquid phase occurs at a crit¬ 
ical magnetic field of 93 kOe.^ However, previous exper¬ 
iments showed that if the magnetic field is applied in the 
easy plane, NTENP underg oes a n Ising quantum phase 
transition at « 113 kOe.l^i^ The high-field phase is 
gapped and antiferromagnetically ordered. Due to its ex¬ 
perimentally accessible critical field and small interchain 
interactions, NTENP in a magnetic field applied in the 
easy plane appears to be a good candidate for the present 
study. 



Figure 1. Crystal structure of NTENP. The magnetic Ni^+ 
ions are surrounded by octahedra with N®“ ions in the equa¬ 
torial vertices. This leads to a single-ion anisotropy with a 
hard axis along a. The chains run along a, where the Ni^^ 
ions are linked by NO^ groups. Perchlorate anions between 
the chains are omitted in the figure for clarity. Two unit cells 
are indicated by dotted lines. 

In this paper, the focus is on the local or q-integrated 
dynamic structure factor, 

/ OO 

(5f(<),5f(0))e“‘dt, (2) 

-OO 

which we measure in the temperature range between 
0.027 K and 12.5 K, for neutron energy transfers Hu up 
to 1.7 meV. These experimental results are used to test 
the predicted scaling law, and to probe the limits of its 
applicability. 


II. EXPERIMENTAL DETAILS 

Half-deuterated (i.e. about 50 % of the hydrogen 
atoms replaced by deuterium) single-crystal samples of 
NTENP, from the same batch as the samples used for the 
experiments of referenc^^^, were characterized by single¬ 
crystal x-ray diffraction. The data were collected on a 
Bruker AXS single-crystal x-ray diffractometer employ¬ 
ing a Cryostream 700 from Oxford Cryosystems for tem¬ 
perature control between 300 K and 100 K. The samples 
were additionally studied in calorimetry experiments. We 
employed a commercial Quantum Design physical prop¬ 
erty measurement system (PPMS) for collecting specific- 
heat data in the range between 300 K and 1.8 K in zero 
magnetic field. Contributions to specific heat from the 
setup and the grease used for mounting the sample were 
subtracted from the data. The temperature rise for each 
data point was 2 % of the sample temperature. 

Neutron scattering experiments were performed on a 
fully-deuterated single crystal of mass 0.5 g which was 
previously used for the experiments reported in refer¬ 
ence^. The crystal was aligned with the neutron beam 
perpendicular to the crystallographic b axis. Neutron 
data were collected on the three-axes spectrometers IN14 
at ILL and FLEXX at HZB,I^ as well as on the time- 
of-flight multi-chopper spectrometer LET at ISIS. A 
150 kOe (IN14), a 148 kOe (FLEXX) and a 140 kOe 
(LET) vertical magnet (H || b) with dilution insert were 
used. The final neutron momentum kf was fixed to 

1.3 A ^ in the experiments on IN14 and FLEXX. On 
LET the incident neutron energy was set to 2.2 meV. 
With these setups, the full width at half maximum 
(FWHM) of the quasi-elastic line of NTENP was mea¬ 
sured to be 0.12 meV (1N14), 0.10 meV (FLEXX) and 
0.05 meV (LET), respectively. For each experiment, the 
critical magnetic field He was determined independently 
by following the field dependence of the magnetic Bragg 
peak at q = (1,0,0) in the ordered phase at H > He- 

III. EXPERIMENTAL RESULTS 

X-ray measurements on half-deuterated single crys¬ 
tals revealed new incommensurate structural satellite 
peaks appearing in the fcLplane for all integer h at low 
temperatures (Fig.j^. The incommensurate peaks per¬ 
sist at least down to 0.5 K as was observed in the 
neutron scattering experiment on LET. They reside at 
5q = ± [0.00(2), 0.21(3),—0.31(2)] around the commen¬ 
surate nuclear Bragg peaks and appear upon cooling at 
about 180 K as is illustrated in Fig.j^a), signaling a pre¬ 
viously unknown structural transition. The occurrence 
of a phase transition is further confirmed by specific-heat 
measurements (Fig.j^b)). The exact values of the tran¬ 
sition temperature and of the satellite-peak intensity is 
not fully reproducible. Moreover, the non-monotonous 
behavior of the satellite-peak intensity is not understood. 
Nevertheless, below 180 K the superstructure is present 
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in all samples studied. Our data is clearly insufficient to 
understand the details of the phase transition and the 
newly established structure. This issue may be subject 
to a further study. However, for the present study the 
precise nature of the phase transition is not important, 
since we only focus on universal magnetic properties at 
the QCP. As will be discussed below, for us, the only rel¬ 
evant consequence of the transition is the residual spin 
gap at He- 



(7,l,-5) + ^ (0,0.21,-0.31) 


Figure 2. X-ray diffraction data measured in NTENP in the 
vicinity of the (7,1, —5) Bragg peak at 300 K and 120 K. Two 
satellite peaks are clearly visible at 120 K. The inset shows 
the 7fcZ-plane at 120 K where red circles mark the principle 
Bragg peaks of the room-temperature crystal structure. 

The main focus of the present work is on the dynamic 
structure factor at the QCP. A crucial point is an ac¬ 
curate determination of the critical field He- ft can be 
determined by observing the magnetic Bragg peaks due 
to AF order stabilized in the high-held phase by weak 
interchain interactions.^ Since there are two spins per 
lattice period along the a axis (see Fig.[^, the AF Bragg 
peaks occur at integer values of h. To determine i2c, 
at the start of each neutron experiment, the intensity 
of the q = (1,0,0) Bragg rehection was followed versus 
magnetic held at a temperature close to the base temper¬ 
ature. Exemplarily, data collected on FLEXX are shown 
in Eig.|^ The critical held was estimated using an em¬ 
pirical power-law ht (solid line in Eig.|^. The actual 
zero-temperature critical held has to be slightly smaller 
than the magnetic held where the magnetic Bragg peak 
appears. Nevertheless, since the system remains in the 
quantum critical regime as long as the temperature is 
large enough compared to the energy gap, small devi¬ 
ations of the used magnetic held from He are accept¬ 
able. The magnetic held chosen for the experiments are 
112 kOe (ELEXX) and 113 kOe (IN14 and LET) respec¬ 
tively. 

On the time-of-hight spectrometer LET the full low- 
energy spectrum could be obtained in the vicinity of the 
AF wave vector. It was collected at 113 kOe at 0.50 K, 



Figure 3. (a) The temperature dependence of the satellite- 
peak intensity measured in NTENP. (b) The specific heat 
features a maximum at about 170 K. To emphasize the max¬ 
imum visually, a second-order polynomial function was sub¬ 
tracted from the total specific heat. The inset illustrates the 
total specific heat (black, solid line) and the fit (red, dotted 
line). 


2.66 K and 5.09 K. Additionally, the spectrum in zero 
held was measured at base temperature (0.5 K) and used 
as the background for the measurements in held. Respec¬ 
tive data are shown in Fig.[^ On the three-axes spec¬ 
trometers IN14 and FLEXX constant-energy scans were 
performed in magnetic held at energies between 0.15 meV 
and 0.40 meV for various temperatures in the range be¬ 
tween 0.027 K and 12.5 K. On ELEXX, the background 
was measured at 0.15 K in zero held and subtracted from 
the raw data (Fig.|^. In contrast, the background for 
the measurements on INI4 was extracted from the data 
obtained in held. Eor this purpose, for each energy a 
polynomial function was htted to the data at all differ¬ 
ent temperatures simultaneously, not taking into account 
the excitation in the window 0.85 < h < 1.15. The data 
was integrated numerically along h in order to obtain 
the g-integrated or local dynamic structure factor S (w). 
Eigure summarizes all data for S (w) collected on the 
three spectrometers for all measured temperatures and 
energies huj <0.9 meV. 
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Figure 4. Peak intensity of the (1,0,0) Bragg reflection at 
40 mK versus magnetic field, as measured on FLEXX / HZB. 
Above a critical magnetic field Hc\t =40 mK the system is 
magnetically ordered and exhibits a magnetic Bragg peak at 
q = (1,0,0). The solid line is an empirical power-law fit to the 
data with fitted exponent 0.5(2) yielding a critical magnetic 
field of H^\ T ^40 mK ~ 112.3(5) kOe > He. 


IV. DISCUSSION 
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A. Residual energy gap at 


The spectrum at 113 kOe measured on LET (Fig.j^b)) 
clearly exhibits an energy gap, A « 0.2 meV, although 
the system is expected to be gapless. A possible explana¬ 
tion of the appearance of the gap is that the applied mag¬ 
netic field deviates from the actual critical held. How¬ 
ever, no closing of the gap was observed in several mea¬ 
surements between 105 kOe and 120 kOe (not shown). 
Therefore, we conclude that the gap does not close at 
any value of applied magnetic held at all. The fact that 
no residual gap has been observed in previous studied^ilSSI 
is most certainly due to an insufficient energy resolution 
in the chosen conhgurations of those experiments per¬ 
formed on three-axes spectrometers. 

The most probable cause for the residual energy gap 
close to the critical held is the observed structural phase 
transition reported in this work. One possible effect is the 
formation of additional anisotropies like Dzyaloshinskii- 
Moriya (DM) interactions, an anisotropic g-tensor or 
in-plane anisotropies evoking a staggered magnetic held 
which prevents the gap from closi ng like in the case of 
the related compound NENP.I^Sllll Indeed, nuclear mag¬ 
netic resonance (NMR) experiments on NTENP revealed 
a held-induced inhomogeneous magnetization for a mag¬ 
netic held applied in the a5*-plane.l^ 


Figure 5. The low-energy spectra of NTENP at 0.5 K mea¬ 
sured on LET / ISIS are shown in zero field (a) and at 113 kOe 
(b). The zero-field spectrum was used as the background for 
the measurements in field. In magnetic held, the S — 1 triplet 
states are split due to the Zeeman effect and the excitation en¬ 
ergy of the lowest excited state is reduced whereas one mode 
persists at Aq ~ 1.2 meV. The residual energy gap which can 
be observed at 113 kOe is discussed in section IIV Al 

B. The scaling law and its empirical range of 
applicability 

Close to a QCP the ^-integrated dynamic structure 
factor should not depend on any properties of the system 
under investigation but is determined by the underlying 
universality class. The only energy scale, for example, is 
provided by the temperature T itself. As a result, it can 
be written in a scaling form: 

5(a;)=T-“$(^), (3) 

where the scaling function <i> and the scaling exponent 
a are specihe to the universality class.^^ This formula 
implies that all ^-integrated neutron data measured at 
different temperatures should collapse on a single curve 
given by d* when plotted versus uj/T and scaled by a 
factor T“. 

In practice, one has to exercise caution. The residual 
energy gap implies that the system is never at the QCP 
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Figure 6. Background-subtracted constant-energy data mea¬ 
sured in NTENP at JY = 112 kOe si on FLEXX / HZB 
at 68 mK around the AF wave vector h = 1. A constant 
offset of 100 counts/~5min is applied between data measured 
at different energy transfers. Solid, black lines are guides to 
the eye. 



/ (kgT) 

Figure 7. Summary of the neutron-scattering data for 
NTENP at the critical magnetic field obtained on IN14, 
FLEXX and LET. The g-integrated dynamic structure fac¬ 
tor is shown versus the temperature-scaled energy transfer 
huj/ (ksT). 


in the experiments of the present work. Despite that, the 
properties of the QCP can still be probed at temperatures 
or energy transfers (frequencies) that are large compared 
to the gap energy such that the system is in the quantum 
critical regime, i.e. 

(keT A) or {huj » A). (4) 

An additional upper limit for temperature and energy is 
due to the presence of the non-critical spin fluctuations 
at higher energy transfers, Aq ~ 1.2 meV, as shown in 
Fig-ib). At too large temperatures and energy trans¬ 
fers, this higher mode is excited and the corresponding 
scattering cannot be separated from the critical part. 
Thus, the second scaling condition to be fulfilled is 

keP,/ku <C Aq. (5) 

Only data fulfilling the scaling conditions, equations 
and (§, can be considered for testing the scaling rela¬ 
tions. The practical upper and lower boundaries were 
determined empirically to be 

(T > 2.9 K) or {Huj > 0.25 meV), (6) 


(T < 7.5 K) and {fiw < 0.9 meV). (7) 

In order to determine the scaling exponent a, the qual¬ 
ity of the overlap between data sets collected at different 
temperatures has to be quantified. For this purpose, an 
approach similar to the one in referenc^^ was used. As 
described in detail in the supplementary material,!^ a x^~ 
like measure V (a) is introduced for the data mismatch 
in the scaling plot produced with a particular scaling ex¬ 
ponent a. Figure]^ shows the scaled data for three differ¬ 
ent scaling exponents a as well as the mismatch function 
V (a) for a in the range between 0 and 1.5. It is minimal 
at ttmin = 0.77(2) implying that the best data overlap is 
achieved for this value of the scaling exponent. 


C. ID Ising QCP 

The observed scaling form for the local spin correla¬ 
tions at the QCP in NTENP can now be compared to 
specific predictions for the ID Ising QCP. For this model, 
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Figure 8. Scaled g-integrated dynamic structure factor S (w) 
measured in NTENP and plotted with scaling exponents 
a = 0, 3 and amin- The scaling exponent with the best data 
overlap, Qmin = 0.77(^ was determined by minimizing the 
function T) (a) (inset) “ For the data with a = 0 an addi¬ 
tional constant prefactor is applied in order to better discern 
the data with 0 = 0 and amin- 


the transverse dynamic structure factor 


^ (g, w) oc T <1 
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( 8 ) 


with the critical exponent aismg = 0.75 and the scaling 
function given by 



The experimentally identified scaling exponent for 
NTENP, a = 0.77(2), agrees remarkably well with this 
expectation. Figure shows the scaling plot for NTENP 
with the Ising scaling exponent oising = 0.75. Addition¬ 
ally, the solid line represents a fit of the scaling func¬ 
tion given in Eq. ( |10[ ) with an overall prefactor (vertical 
shift in the log-log plot) as the only adjustable param¬ 
eter. Once again, an excellent quantitative agreement 
with the data is apparent. 



hco/ i\T) 


where q is the momentum transfer relative to the critical 
wave vector, u is the spin wave velocity, P is the gamma 
function and 5 denotes the imaginary part of the argu¬ 
ment. Integrating over q yields the scaling form. 


Figure 9. The scaling plot for the local dynamic structure 
factor measured in NTENP with the Ising scaling exponent 
aising = 0.75 is depicted. The black line corresponds to the 
predicted universal scaling function for the quantum critical 
point of an Ising chain in a transverse magnetic field as de¬ 
scribed in the text, Eq. (101. 


^(cc) = r-tci>(^). 


(9) 
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V. CONCLUSIONS 

In summary, the scaling relation for the local dynamic 
structure factor of a linear-chain antiferromagnet has 
been verified at a field-induced Ising QCP over three 
orders of magnitude in huj/ {k^T). The use of neutron 
spectroscopy and the coverage of a wide range of energy 
transfers allows to evade non-critical spin fluctuations. 
It also provides a way to avoid the unwanted terms in 
the Hamiltonian that prevent an exact realization of the 
ID Ising QCP in any realistic material. In this respect, 
neutron spectroscopy shows more flexibility than other 
techniques, e.g. NMR, which essentially measures the 
same local dynamic structure factor but is restricted to 
a single, low measurement frequency.^^^^ 
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I. QUANTIFYING THE DATA OVERLAP 

For the quantification of the quality of the data overlap for different scaling exponents a the following parameter 
was used 


D(a) = 




Aa)-a„ 




(a))" 


^n.m ^in 




( 1 ) 


where n and m run over all different experiments and temperatures, and jm run over the energies for experiment 
/ temperature n and m, respectively, an,i„ (a) = TffS{uJn,i„,Tn), and Aan,i„ (a) is the corresponding error. The 
prefactor selects data in the proper uj/T range. 
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